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Abstract 

It is shown that, by appropriately defining the eigenfunctions of a function defined on the 
extended phase space, the Liouville theorem on solutions of the Hamilton-Jacobi equation can 
be formulated as the problem of finding common eigenfunctions of n constants of motion in 
involution, where n is the number of degrees of freedom of the system. 


1 Introduction 

In the framework of the Hamiltonian formulation of classical mechanics, the Liouville theorem asserts 
that, for a mechanical system with n degrees of freedom, if we have n constants of motion in 
involution, Fi, F2, ..., Fn (that is, {Fi,Fj} = 0 for i,j = where { , } is the Poisson 

bracket), then a complete solution of the Hamilton-Jacobi (HJ) equation can be found by quadrature 

More precisely, ii Fi{qi,pi,t), ... ^ Fn{qi,pi, t) are n constants of motion in involution, that is, 

QF 

^ + {C,iJ} = 0, z = l,2,...,n (1) 

and 

{Ft,Fj} = 0, i,j = 1,2,...,n, (2) 

where { , } denotes the Poisson bracket (with the convention {qi,Pj} = Sij), then, assuming that 

so that, locally at least, we can express the pi as functions of qj,Fj, and t, the differential form 

Pi (Qj ,Fj,t)dqi- H{qi,p^{qj,Fj,t),t)dt (4) 

is the differential of some function, S{qi,t), which is a complete solution of the HJ equation (here 
and in what follows, there is summation over repeated indices). 

The aim of this paper is to show that the Liouville theorem can be formulated in another form, 
closer to the standard formalism of quantum mechanics. Specifically, we shall show that if S{qi,t) 
is a common eigenfunction of the functions Fi, F 2 ,..., Fn (a concept to be defined below) then, by 
adding to S an appropriate function of t only, one obtains a complete solution of the HJ equation. 

In Section 2 we present the definition of the eigenfunctions of a function f{qi,Pi,t) and we 
show that two functions, f{qi,Pi,t) and g{qi,pi,t), have common eigenfunctions if and only if their 
Poisson bracket vanishes. Then, we prove that, if conditions (HJ-® hold, a common eigenfunction 
oi Fi, F 2 ,..., Fn is, up to an additive function of t only, a complete solution of the HJ equation. In 
Section 3 we give some illustrative examples, emphasizing the fact that we can make use of constants 
of motion that depend explicitly on the time. 

The statement of the Liouville theorem presented here allows us to see that the Liouville theo¬ 
rem is analogous to one of the methods employed to solve the Schrodinger equation, where we look 
for the common eigenfunctions of a complete set of mutually commuting operators that also com¬ 
mute with the Hamiltonian (e.g., for a spherically symmetric Hamiltonian we consider the common 
eigenfunctions of and H). 
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2 Eigenfunctions of a function and complete solutions of the 
Hamilton—Jacobi equation 

We start by giving the definition of the eigenfunctions of a real-valued function f{qi,Pi, t): We shall 
say that S{qi,t) is an eigenfunction of f{qi,Pi, t), with eigenvalue A, if S' is a solution of the first-order 
partial differential equation 

aq 

/(..^,t) = A. (5) 

(It may be noticed that if / is a time-independent Hamiltonian, then ([5]) is the corresponding time- 
independent HJ equation.) We note that if S{qi,t) is an eigenfunction of f{qi,Pi,t) with eigenvalue 
A, then so it is S{qi, + for any function (j){t) of t only, and that the solutions of ([5]) will depend 
parametrically on A. Of course, in order for ([5]) to be a differential equation, / must depend on one 
of the Pi, at least. 

For instance, according to this definition, the eigenfunctions of the function 

F{q,p,t) = mujqsmujt + pcosujt, (6) 

where m and uj are constants, are the solutions of the differential equation 

• .93 ^ , 

mujqsmujt + cosujt = A, 
oq 


which can be readily integrated giving 

S = Agsecwt — tanwt-I-(7) 

where 4>(t) is an arbitrary function of t only. Note that A may be a function of t. 

If S{qi,t) is a common eigenfunction of f{qi,Pi,t) and g{qi,pi,t), with eigenvalues A and p, 
respectively, that is, S satisfies and 


. dS , 

9{qt, = P, 


then, differentiating with respect to making use of the chain rule, we obtain 


hence, 


df df d^S ^ , dg dg d^S 

a —^— a —® 'a—^— a —® 

^Pj ^Pj 


if, 9 } 


df dg 

dg df 

df d^S dg 

dqi dpi 

dqi dpi 

dpj dqjdqi dpi 

/ 925 

d^S \ 

df dg 

\dqjdqi 

dqidqj J 

dpi dpj 


+ 


dg d^S df 
dpj dqjdqi dpi 


Thus, if f{qi,pi,t) and g{qi,pi,t) possess common eigenfunctions, then {f,g} = 0. 

In order to see that the converse is also true, we now assume that {/, 5 } = 0. If / and g are 
functionally independent, then there exists, locally at least, a set of canonical coordinates, Qi,Pi, 
such that. Pi = f and P 2 = g. Then, the eigenvalue equations for / and g are dS/dQi = A and 
dS/dQ 2 = p, which have the simultaneous solutions S = AQi -I- pQ 2 + , Qn, t), where (p is 

an arbitrary function of n — 1 variables, thus showing that / and g have common eigenfunctions. 
(Note that this does not mean that every eigenfunction of / is an eigenfunction of g (cf. [7], Sec. 
2.9).) The expression for S in terms of the original coordinates, {qi,t), is not given by the simple 
substitution of the Qi as functions of {qi,pi,t) [ 8 ]; what is relevant here is the existence of common 
eigenfunctions for / and g. 

In the case where / and g are functionally dependent, the eigenvalue equations for / and g are 
equivalent to each other and, trivially, possess common solutions. 
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2.1 Alternative formulation of the Liouville theorem 

We now assume that Fi,..., Fn are n functions satisfying 0-0, and we consider a common 
eigenfunction S{qi,t) oi Fi,..., Fn, with eigenvalues Ai,..., An, respectively, then, assuming that 
the eigenvalues are constant, differentiating with respect to t both sides of the equation 

ae 

F,iq„—,t) = X„ ( 8 ) 

making use of the chain rule, the Hamilton equations and o, we have 


0 = 


dFi . OF 

% ^ 
c 

dqj dpj c 
dF, { d^S 


dS 

dF, 






:dqj 

dt 






/ d^S 


d^S ^ \ 

dF, 

dH 

, dF, 

dH 

V dtdq- 

i ' 

dqkdqj^*") 

1 ^ 

dpj 

h 

dqj 

d^S 

dH 

dH\ 



. .,n. 


dqkdqj 

dpk 

1 

1- 

i — 

: i, . 



By virtue of @, the last equations are equivalent to 


0 = 


d‘^s 

dtdqj 

A _ 

dqj [ dt 


+ 


dS 


d^S dH dH 

dqkdqj dpk dqj 

rr, dS ' 

H{qk,j—,t) 

dqk 




which implies that the expression inside the brackets is a function of t only. 


dS dS , 


Thus, 


S = S- 


X{u )du 


is a solution of the HJ equation. We can verify that this solution is complete by differentiating (|8]) 
with respect to Aj, which gives 

dF, d^S 
dpkdXjdqk 

Taking into account dS]), this last equation shows that dei{d^S/dXjdqk) ^ 0. 


3 Examples 

In this section we give some examples of the method presented above. 


3.1 One-dimensional harmonic oscillator 

The function 

F(( 7 , p, t) = rtiLoq sin ut + p cos ut 

already considered above [see ®], is a constant of motion if the Hamiltonian is given by 

where w is a constant. According to the results of the preceding section, if A is a constant, 

mui 2 , / N 

b = Xq sec ujt - —q ta.n ojt + (p{t) 


(9) 


( 10 ) 
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must be a solution of the HJ equation for the Hamiltonian (jH]), if the function ^ is appropriately 
chosen. A direct computation yields 


1 


dS 


mUp 2 2 /// \ 


2 m \ dq 

and, therefore, choosing (j){t) = —tanwt/2ma;, we obtain the complete solution of the HJ equation 


, muj^ r 

b = Xqsecujt — - -1- —q" 

2 m 2 


tanwt 


w 


Note that, in this case, H is also a constant of motion and, as pointed out above, the equation 
that determines the eigenfunctions of H is just the time-independent HJ equation. However, the 
constant of motion ([0]) leads to simpler expressions. 


3.2 Particle in a time-dependent force field 

As a second example we consider the time-dependent Hamiltonian 

H = - - ktq, 

Zm 

where fc is a constant. One can readily verify that 

kt^ 

F = P- — 

is a constant of motion and that the eigenfunctions of F, i.e., the solutions of 

dS ke 


dq 


= A, 


are given by 


kt^ 

S = Xq+ —q + (l>{t), 


where (j){t) is an arbitrary function of t only. Then 


1 


2 m \ dq 


dS 


2 dt 2m 


XkF 

2m 


2+4 


-b 


kH 

8 m 




hence, choosing 


,, . A^t Xkt^ k'^t^ 
- 


2m 6m 40m ’ 
m is a complete solution of the HJ equation (which is not separable). 


( 11 ) 


3.3 Particle in two dimensions 

As a final example we consider the Hamiltonian 


H = 


1 

2 m 


Px 


eB 


y] +\Pv- 


eB 


( 12 ) 


which corresponds to a charged particle of mass m and electric charge e in a uniform magnetic field 
B. The functions 


„ 1/1 \ I ■ . muj , . 

Pi = — [I + cosojtjpj; ——py sinu!t——xsmojt - —{l — cosujt)y, 

„ 1/1 \ I • /I 1 

p 2 = — [1 + COSojtjpy +—px Sinu!t-\ — —ysinojt-l —^(1 —cosa;t)x. 


(13) 

(14) 


where w = eBjmc^ are constants of motion in involution, which correspond to the values of the 
canonical momenta p^ and Py, respectively, at < = 0. 
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From (I13|) and (fHl) one finds that the common eigenfunctions of Fi and F 2 , with eigenvalues Ai 
and A 2 , respectively, are 


S'= Aicc + A 22 /+ tan X 2 X — Xiy - ^{x'^ + y^) 

where is an arbitrary function of t only. Substituting this expression into the HJ equation one 
finds that S is a solution of this equation if and only if 


Ai^ + A2^ 
2m 


sec^ ^ujt + = 0, 


hence, 

S = Aia: + X2y — 

is a complete solution of the HJ equation. 


/, mw \ 2 / mw \ 

(Ai + —y) +(A2-—x) 


tan 


4 Concluding remarks 

As pointed out above, the formulation of the Liouville theorem given here makes use of terms anal¬ 
ogous to those employed in the standard formalism of quantum mechanics, thus providing another 
example of the parallelism between both theories. Another advantage of the version of the Liouville 
theorem given above is that its proof is shorter than those usually presented in the textbooks. 
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